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Nowadays, direct torque control (DTC) of AC drives is one of the most 
interesting fast dynamic response motor control techniques. Permanent 
magnet synchronous motors (PMSM) are now replacing induction motors 
(IM) in a wide range of applications because of their salient merits such as 
high efficiency and compactness. Combining DTC and PMSM in one drive 
may result in excellent performance not only for torque dynamics but also for 
efficiency. DTC relies on the estimation of the stator flux and torque to be 
compared with the set points of the torque and flux. In the literature, there are 
various flux estimation methods and each one has its own merits and demerits. 
Comparisons of these various methods are discussed in this paper showing 
their effect on the drive performance through MATLAB/Simulink 
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1. INTRODUCTION 

Direct torque control (DTC) of AC drives may be considered as one of the most interesting fast 
performance control methods nowadays. Permanent magnet synchronous motors (PMSM) are now replacing 
induction motors (IM) in many applications because of their valuable characteristics especially compactness 
and high efficiency. The inclusion of DTC and PMSM in one drive system may result in high dynamic 
response [1]. 

DTC offers high dynamic torque and flux with simple controller that uses two hysteresis controllers 
and a switching table to select the most suitable voltage vector to meet the torque and flux requirements. No 
machine parameter dependences, no voltage modulation techniques, no axes transformation needed, and no 
current controllers are the advantages of DTC over FOC. Although, DTC gives excellent dynamic responses 
for both torque and flux, the ripple content is high [2], [3]. 

For high performance drives, field-oriented control (FOC) or direct torque control (DTC) may be 
utilized. Based on [4] and [5], comparing FOC and DTC it is concluded that DTC provides a simple control 
drive with a better dynamic response than FOC as in DTC, no motor parameters dependence, no position 
sensors, no current control loops, no need for reference frame transformation, and voltage modulation is not 
required. DTC has hysteresis controllers with switching look up table. DTC gives separated control for torque 
as well as flux. It chooses the most suitable inverter voltage vector from a pre-designed look up table which 
would result in fast torque response. The selection is conducted to limit the flux and torque errors within their 
hysteresis bands. The DTC has the demerits of high torque and flux ripples, variable inverter switching 
frequency, and bad performance at low speeds [4], [5]. Moreover, many researches have been conducted to 
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solve DTC problems. Some studies have been made that neither use switching table nor hysteresis controllers 
but they calculate the voltage vector and realize it by one of the voltage modulation methods like Space vector 
modulation (SVM-DTC) [6], [7] or Sinusoidal pulse width modulation (SPWM-DTC) [8]. Whereas in [9], 
[10], the voltage vector is applied to a certain ratio of the sampling time, which reduces the torque ripples 
remarkably in a control method called duty ratio control. Model predictive control (MPDTC) also replaces the 
table and the hysteresis controllers by a prediction algorithm and a cost function. The prediction algorithm uses 
the motor equations to find the trajectory of the torque in the next sampling times that also decreases the torque 
ripples [11]. Okumuş and Aktaş [12] presented the hysteresis band of the hysteresis controllers is adapted in 
real time with the stator flux and torque errors variation. 

Moreover Tarmizi et al. [13] presented the DTC of Brushless DC motor is proposed by injecting the 
triangular waveform and using PI controller to reduce the torque ripples and obtain constant switching 
frequency. Recently, nine-switch inverter is utilized in [14] to drive two independent induction motor loads 
with DTC to reduce ripples in torque and flux with higher efficiency. While Echeikh et al. [15] presented five- 
phase induction motor drive using backstepping control with constant switching frequency is developed. 
Whereas, improvement of DTC applied to doubly fed induction motor under variable speed is proposed in [16] 
using multilevel inverter. Furthermore, SVM-DTC strategy of PMSM employing ultra-sparse matrix converter 
that uses only nine switches is developed in [17] to reduce the flux and torque ripples with a cheaper converter. 
While Nasr et al. [18], torque ripple suppression for interior permanent magnet synchronous motor IPMSM 
using finite element analysis based model predictive-DTC is derived to optimize the duty ratio simultaneously 
with the voltage vector selection to guarantee further torque ripple reduction under steady-state operation, 
especially at low speeds. Although these relatively new methods for DTC can solve various problems of the 
conventional method; these new techniques require the use of position sensors which increases system 
complexity, cost, and reduces reliability. Machine parameters dependence is also a disadvantage of these new 
techniques, thus on-line parameter estimation, temperature rise, and saturation should be considered. Increased 
computations and control algorithm complexity need the utilization of powerful processors that have high 
performance capabilities, therefore cost is affected [5]. The conventional DTC drive controller basically 
consists of flux estimator to deduce the torque. Any error in flux vector estimation rather than its real value 
would result in errors in both flux and torque control loops. Therefore, flux estimation is a crucial stage in a 
DTC drive [19]. 

In this paper, four different methods of flux estimation are compared through MATLAB/Simulink 
simulations to interpret their effect on the DTC drive performance. The paper is organized as following; 
section 2 provides the modeling of the PMSM motor in the stationary reference frame. Section 3 discusses the 
concept of DTC as applied to PMSM. While in section 4, the control strategy of the flux vector is presented. In 
section 5, the voltage selection looks up table is outlined. The various flux estimation methods are discussed in 
section 6 while the results and discussions are provided in section 7. Finally, conclusions are given in section 8. 


2. MODELING OF PMSM IN THE STATIONARY REFERENCE FRAME (D-Q PLANE) 
The stator voltage space vector (us) of PMSM in the stationary reference frame can be expressed by (1). 


dws 


us = Rsis + ae (1) 
Where: (Rs): the stator resistance and (¥,): the stator flux linkage. The stator flux linkage can be given by (2). 
ws = Lsis + YPM e/) (2) 


Where: (Ls): the stator self-inductance, (¥pm): the permanent magnet stator flux linkage, and (©,): the rotor 
angle. Substitution (2) into (1) results in: 
di 


Us = Rs is + Ls Z$ +j Or Yem e) °r (3) 


de 
Where œ, = ye 


The (3) represents the motor voltage equation in the stationary reference frame. This equation can be 
decomposed into its D-Q components as in (4) and (5): 


while Ypy and L, are constants. 


usD = RsisD + LsD “ae — wr YPM sin(Or) (4) 
usQ = RsisQ + LsQ Ssg + wr YPM cos(Or) (5) 
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3. CONCEPT OF DTC AS APPLIED TO PMSM 
It is known that quick torque response is obtained by changing the position of the stator flux linkage 
space vector with respect to the rotor flux space vector as (6) shows. 


Te ==P Ce) [vs]. pm sin (5) (6) 


The stator flux is maintained constant within certain limits to ensure the operation at the rated condition and to 
avoid the saturation region. Therefore, the control of torque is obtained by controlling the torque angle (6) as 
the permanent flux (Wp) is also constant. 

In the stationary reference frame, the stator flux linkage space vector is the integration of the stator 
EMF and itis obtained by in (1). If the stator resistance voltage drop is neglected for simplicity, then the stator 
flux is the integration of the applied voltage. Thus, in a small time period, the increment of the stator flux relies 
on the applied voltage. Then, the desired trajectory of the stator flux can be achieved by selecting the 
appropriate inverter voltage vector. The following analysis shows this idea [20]: 


Š dy 
Us = Rg is + oT then, 


dYs =u- R; is then, 
dt 
ws = f (us — Rs is) dt 7 


It is shown from (7) that the stator flux linkage space vector ( w, ) can be controlled by controlling the stator 
voltage ( us ). Where, ( us ) is one of the VSI eight voltage vectors. For simplification, the resistive voltage 


drops (Rs is) can be neglected, then (8). 
y, = fu, dt (8) 
Also, for a small period of time (At), the voltage vector can be assumed constant, then (9). 


(9) 


gn ee 


t=0 


Where, Ys : is the initial value of the stator flux linkage. For a switch-on time (At), the increment of the 
—lt=0 


stator flux can be obtained by (10). 


Aw, =u, . At (10) 


This means that the stator flux linkage space vector moves by Aw, in the direction of the applied stator voltage 
vector, as shown in Figure 1. 


Us : us 


Figure 1. Inverter voltage vectors and the corresponding flux variation in time At 


4. CONTROL OF STATOR FLUX-LINKAGE SPACE VECTOR MAGNITUDE 

To choose the appropriate voltage vectors required to control the magnitude of the stator flux space 
vector, the voltage space vector plane (D-Q plane) or (a-f plane) is divided into six sectors K(1) through K(6), 
Figure 1, such that, [20]: 
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(2N - 3) È < K (N) < (2N-1)? (11) 


where, N = 1, 2,3...... , 6. The idea is to make the switching frequency as low as possible so that the most 
suitable voltage vector is the one requires the minimum number of switching and drives both the stator flux 
and the torque errors in the needed direction. In the k" sector, the stator flux can be increased by applying one 
of the: k, k+1 (anticlockwise), or k-1 (clockwise) vectors. Therefore, flux can be increased by applying the 
voltage vector belongs to the sector or any of the adjacent voltage vectors. To decrease the stator flux, one of 
the voltage vectors is applied: k+2 (anticlockwise), k-2 (clockwise), or k+3. A deviation by a hysteresis band 
equal to tA¥, is allowed for the actual flux with respect to the reference flux [20]. The block diagram of the 
DTC scheme is given in Figure 2. 


3- Level Hysteresis 
Comparator Voc 


- i 
3Phases |e 
a ee 
imati. 2Phases 
Estimation 
Voa Vasc 


Figure 2. The block diagram of DTC of PMSM 


5. THE VOLTAGE SWITCHING LOOK-UP TABLE 

It is a pre-designed table as shown in Table 1 that utilizes its inputs from the torque and flux 
comparators and uses the sector number to achieve the needed voltage vector. Its criterion is based on the 
previous idea of the stator flux vector control [19]. Assume that the flux vector is in sector (1) and dw, = 1& 
dr =1, then both of torque and flux are needed to be increased. Thus, uz (110) is applied, Figure 3. Note that u: 
(100) cannot be applied as torque will be reduced as the torque angle will be decreased. Now, dw, = 1& dr =0, 
the flux does not be changed instantaneously and torque has the priority to be changed, then a zero-voltage 
vector is applied that gives minimum switching states, u7 (111). Now, dw, = 1& dr = -1, flux is required to be 
increased and torque is required to be decreased, then us (101) is applied. In the case of dw, = 0 & dr =1, u3 
(010) is applied to increase the torque [20]. 


Table 1. The voltage switching B E- T+ FL T+ 
look-up table uz m2. 
d, d KU) K() KG) KA KG) K(@) Pa 
1 1 Uy u u4 us U6 u “ue a“ sector 
u7 Uo u7 Uo u7 Uo 
-1 U6 uy U2 U3 U4 Us 
0 1 us uy us U6 u Ww uš sc NG 
uo u uo u uo u F- T- FRE 
-1 Us U6 uy Uo U3 U4 


Figure 3. Voltage vector selection concept 


6. VARIOUS STATOR FLUX ESTIMATION METHODS 

Comparisons between some implementation methods of the flux estimation are conducted. The 
simulation results will help us know the advantages and disadvantages of each method. Through our discussion, 
problems encountered in each method will be stated and suggested solutions against them will be given too. 


6.1. Pure integrator flux estimator 

This method is a direct application of in (7). This means an open loop integration of the induced EMF 
(Vs —is Rs) so, voltage and current measurements are needed. It is simple and does not rely on motor parameters 
except the stator resistance. Also, it does not need the rotor position. However, the implementation of the pure 
integrator has some problems as: DC drift problem, stator resistance variation, initial value problem, operation 
at low frequencies. 
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6.1.1. DC drift problem 

Errors in voltage or current measurements lead to a drift in the estimated flux value. Although it may 
be small, the integrator may finally get into saturation. These errors may be resulted from the following: i) 
Sensors may result in phase shift, ii) Added gains or scaling factors, iii) Quantization errors in digital systems. 
Figure 4(a) indicates estimated stator flux linkages when dc offsets are present. Speed ripples are shown in 
Figure 4(b) when offset exists. Figure 4(c) and Figure 4(d), show flux and speed without any offsets in 
measurements [21]. 


Speed, rpm 
z 


Time, see. 


(d) 


Figure 4. The error effect [21] (a), (b) with offset on flux and speed and (c), (d) without offset on flux 
and speed t 


6.1.2. Stator resistance variation 

An accurate stator resistance value has to be utilized to ensure the correctness of the estimated stator 
flux. The stator resistance changes with the temperature variation and the frequency (skin effect), so it must be 
adapted to these changes. Stability may be affected if the stator resistance value that is used in the control 
algorithm is different from the actual resistance. Figure. 5 represents some results of flux response when the 
stator resistance is increased by 30%, and the controller uses the nominal value [21]. Online estimation of the 
stator resistance or using of thermal machine model is an eligible solution to this problem. A stator resistance 
estimator was proposed in [21] utilizing the relation between current and resistance variations. 


6.1.3. Initial value problem 


Unlike IMs, the initial value of the stator flux vector w, 


is not zero in PMSMs due to the existence 
0 


of the PM but it can be obtained using the rotor position. Therefore, the rotor position should be measured 
using position sensor or estimated using a sensor less technique to ensure correct startup of the drive. If the 
drive is started with a false initial stator flux value, a drift may be presented in the flux different from its 
accurate trajectory [22]. 


6.1.4. Operation at low frequencies 

The pure integrator becomes inaccurate at low speeds. This is because the stator voltage becomes very 
small and both the resistive voltage drop and the inverter switches voltage drop cannot be neglected. The pure 
integrator performs well at high frequencies (speeds) since the resistive voltage drop is very small and the 
voltage drop across the inverter switches may be ignored [20]. 


6.2. Low pass filter (LPF) flux estimator 


A low pass filter (LPF) can be utilized instead of the pure integrator to solve the problems of the drift 
and the initial value problem. Starting from the motor voltage in (12). 
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vs = Rsis + pws (12) 
where p = d/dt, "the derivative operator”. 
A low pass filter is obtained by replacing: 
1 T 
p > 1+pT (13) 
(14) 


4 1+pT 
And hence,vs = Rs is + z5 ws 


where T is the filter time constant and the cutoff frequency œc = : . The cutoff frequency should be selected 
carefully. It should be small enough to act as a pure integrator for a wide speed range. The transfer function of the 


low pass filter is H(s) = —— = — , Which means that a damping effect by a time constant of "T" will help in 


solving the initial value problem but a small steady state error will be introduced in the stator flux final value [23]. 


YQ, Wb 


sQ, Wb 


“08 -04 o4 o6 


02 O 02 
Wp, Wb 


(b) 


Figure 5. Flux trajectories with (a) correct resistance value and (b) 30% increased resistance value [21] 


6.3. Modified integrator flux estimator 
Hu and Wu [24], three modified integrators were proposed to estimate the stator flux over a wide 
speed range. The output of these integrators can be represented by (15). 
wW 
yu a (15) 


St+We St+We 


Where “x” is the integrator input and “z” is a compensation signal. The modified integrator acts as a LPF if the 
compensation signal equals to zero. If the compensation is taken from the integrator's output, the modified 
integrator acts as a pure integrator. The three integrators presented in [24] differ in the way of the compensation 
signal. From the three methods, the third one provides the best performance. It will be discussed in some detail. 

The third algorithm is designed for motor drives that operate at variable flux. The algorithm relies on 
the fact that the flux vector should be orthogonal to the EMF (which is the stator voltage minus the resistive 
voltage drop). The orthogonality is tested by obtaining the dot product of the back EMF and the flux linkage. 
The compensation signal amplitude is controlled using a PI controller. The algorithm block diagram is 
displayed in Figure 6 [24]. 

The estimated flux vector “w,” is the sum of two vectors; a feed-forward vector "w, " and a feedback 


vector "wz" . If there is neither drift nor initial value problems, the orthogonality will be achieved and the dot 


product of the two vectors will equal to zero. If a drift or initial value problem exists, the flux and EMF will 
not be orthogonal, and the angle "y" becomes greater than 90°. This leads to an error signal Ae that is (16). 


hea SEEM a SEMEL (en RORY (16) 


ws] Isl 


Where " y" is the angle between the flux and EMF. 
The output of the PI controller can be expressed as (17). 


Yomp = (kp + “£).A (17) 
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Where kp and kj are the PI controller gains. 
For example, if the value of "», " is increased to " w, \" as given in Figure 7 due to an offset or initial 


value problem, the angle "y" will be greater than 90°. A negative error will be introduced. The PI controller 
output "Yemp " and the feedback vector will be decreased. Therefore, the flux vector "Ys \" changes its position 


toward the "original 90°-position" until the orthogonality is occurred. The opposite case will happen if "y 
becomes less than 90°. The transient as well as the steady state responses of this modified integrator are 
relatively good. It has a problem that is it needs carful tuning for both the PI controller and the low pass filters 


so that complexity is increased. 
1 Pp: Ysp 
T sta, ? 
EMFp 


EMF 


Figure 6. Modified integrator with Figure 7. Vector diagram of the 
adaptive compensation modified integrator 


6.4. Tunable low pass filter flux estimator 

Singh et al. [25], a low pass filter with tunable cutoff frequency is presented for flux estimation for 
IM drive and it offered good results. The flux distortion and ripples are decreased, no flux drop at sector 
boundaries, and the flux trajectory is a pure circle which implies the effectiveness of this estimator over the 
conventional constant cutoff frequency LPF mentioned in section (VI.B) above. As mentioned before, a LPF 
is an eligible solution to the problems of drift and initial value problem that may occur when utilizing a pure 
integrator, but flux distortion may be introduced when the cutoff frequency is higher than the synchronous 
frequency. Thus, the cutoff frequency value should be carefully chosen [25]. A low pass filter with a tunable 
cutoff frequency relying on the synchronous frequency can be a good solution. It is depicted in Figure 8. 
Utilizing a simple relation between the cutoff frequency œ. and the synchronous frequency œe that is: œ. = K * 
@e; the cutoff frequency can be obtained and K is a constant where: 0.1 < K <0.5. The synchronous frequency 
can be obtained as (18). 

we = Sos ese (18) 


1 
S+ Wc 
Wps-Egs — Eps-¥as 
=n] | « | Pps? + Yos? 


S+ We 


Figure 8. Tunable-low pass filter flux estimator 


7. RESULTS AND DISCUSSIONS 

The DTC drive system will be tested via simulations using MATLAB/Simulink for each flux 
estimation method mentioned above. The used PMSM has 6.1 N.M full load torque and full load speed of 180 
rad/sec electric. Motor parameters are shown in Table 2. The Simulink model is shown in Figure 9. 
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Table 2. Motor parameters 


Parameter Value Parameter Value 
Bus Voltage "Vpc" 300V Stator Resistance "Rs" 14Q 
Pair Poles "P" 3 D-axis Inductance "Ld" 0.0056 H 
Permanent Magnet Flux " py" 0.1546 Wb Q-axis Inductance "Lq" 0.0056 H 
Moment of Inertia "J" 0.00176 Kg.m’ _ Friction Coefficient "B" __ 0.0038 N.m.s 
ei T Error —>|n1 Outi Va — vsph > 
Speed 
T torque controller Vb 
Ve VsQ 
Flux Error >) 
pj ere flux controller abe to DQ 
flux ref gen switching table eae 
Fl (oom) 3 N — flux angle}#— 
| stator flux sector calculation pa Int ia >| 
= a 1387 he Fh flux angle In2 ib >» 
flux calculation = 
o 
ic +| 
Lic flux &torque estimator Dat ab ie 
torque o abc 
AW DQ flux 
lä- Speed Ref 
m bl Torque Ref speed Het icles 
speed + 


Saturation 


pi Controller 


Figure 9. Simulink model for DTC of PMSM 


7.1. Pure integrator flux estimator 

As mentioned before, the flux estimation using the pure integration technique is the simplest method 
but it has some problems as DC drift, initial value, and stator resistance variation problems. In simulation, the 
initial rotor position can be adjusted to be equal to zero; therefore, there are not any measurement errors or 
resistance variation in simulation. Therefore, this method may be considered as the nearest one to the actual 
flux and it will be used as a reference for comparing the other flux estimation methods. The motor is operated 
at rated torque and speed 180 rad/sec in the speed control mode. Stator Flux linkage magnitude, phase, D-axis 
flux, and Q-axis flux, are shown in Figures 10, 11, 12, and 13, respectively. 


7.2. Low pass filter flux estimator 

LPF is expected to solve the problems of the pure integrator but phase and magnitude errors are 
occurred especially at low frequencies. The dynamic performance of the LPF flux estimator is poor but it is 
acceptable to some extent at steady state. Figures 14, 15, 16, and 17 show: stator flux magnitude of LPF 
estimator, flux angle of LPF estimator, D-and Q- axes fluxes of LPF estimator compared to the pure integrator 
estimator respectively. 

It is clear from the above figures that the LPF flux estimator provides poor results at transient state 
but it is acceptable to some extent at steady state conditions. The steady state error of the magnitude and phase 
of the flux is relatively high for the LPF estimator. The selection of the filter cut-off frequency is critical in 
implementing such flux estimator type. 


Stator Flux Linkage Magnitude Stator Flux Angle 
0.23 f iia T T T T T T T T j 4 T T T i ig T 
0.22 | 3 
207 | 2b = ; - aa 
ao} 
$ 0.2 g ik 
20.19 $ o | 
o < 
= 0.18} \ 4 %41 
= ira 
i 0.17; —4 t -2 | | 
ee 3 
0.154 i a L L L = 4h i i L A 
0.02 0.04 0.06 0.08 0.1 0.12 014 016 0.18 0.2 0 0.02 0.04 0.06 0.08 0.1 0.12 014 0.16 018 0.2 
time "sec" Time "sec" 


Figure 10. Stator flux magnitude using pure integrator Figure 11. Stator flux angle using pure integrator 
estimator estimator 
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Figure 12. D-axis flux using pure integrator estimator 
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Figure 14. Stator flux magnitude of LPF compared 
to pure integrator estimator 
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Figure 16. Stator D-axis flux of LPF compared to 
pure integrator estimator 


7.3. Modified integrator flux estimator 
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Figure 13. Q-axis flux using pure integrator estimator 
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Figure 15. Stator flux angle of LPF compared to 


pure integrator estimator 
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Figure 17. Stator Q-axis flux of LPF compared to 
pure integrator estimator 


Figures 18, 19, 20, and 21 present the flux magnitude, flux angle, D-Q axes fluxes for the modified 
integrator flux estimation method compared to the pure integration flux estimation method respectively. It is 
clear from the above figures that the transient response is improved due to the compensation effects of the 
modified integrator. Also, the steady state response of the modified integrator is acceptable and better than the 
LPF estimator. The steady state errors for both flux magnitude and phase are decreased compared to that of the 
LPF estimator. As mentioned before, complexity is the major disadvantage of this method of flux estimation 


and it needs many tunings. 
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Figure 18. Stator flux magnitude of modified integrator Figure 19. Stator flux angle of modified integrator 
compared to pure integrator estimator compared to pure integrator estimator 
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Figure 20. Stator D-axis flux of modified integrator Figure 21. Stator Q-axis flux of modified integrator 
compared to pure integrator estimator compared to pure integrator estimator 


7.4. Tunable low pass filter flux estimator 

Figures 22, 23, 24, and 25 display the flux magnitude, flux angle, and D-Q fluxes of the tunable LPF 
estimator compared to the pure integrator estimator method respectively. The tunable LPF estimator provides 
the best results of the above flux estimation methods in both transient and steady state conditions due to its 
brilliant criteria for tuning the cut off frequency of the filter as mentioned before. Its simple algorithm structure 
is another advantage of this method of flux estimation. 
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Figure 22. Stator flux magnitude of tunable LPF Figure 23. Stator flux angle of tunable LPF 
compared to pure integrator estimator compared to pure integrator estimator 
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Figure 24. Stator D-axis flux of tunable LPF Figure 25. Stator Q-axis flux of tunable LPF 
compared to pure integrator estimator compared to pure integrator estimator 


8. CONCLUSION 

Various methods are proposed in the literature to estimate the flux and the most powerful methods are 
compared in this paper via MATLAB/Simulink simulations. Pure integration method of flux estimation is the 
basic method that has a lot of disadvantages like offset, resistance variation dependence, initial value problem, 
and low performance at low speeds. Low pass filter (LPF) is designed to solve the pure integration problems 
but its performance at transient state is poor and there is some drift in the flux from the actual flux. The modified 
integrator is also one of the interesting flux estimation methods that presents better performance at both 
transient and steady state conditions than that of LPF estimator but it contains PI controller and low pass filters 
so many tunings are needed and complexity is increased. A flux estimator uses tunable-cutoff-frequency-LPF 
is simple and accurate. It is utilized for IM before but never be used for PMSM. It is also compared and gives 
excellent results in both transient and steady state periods which implies the effectiveness of it for PMSM as 
well as IM. Therefore, in this paper, the different flux estimators were compared via simulations and the 
conclusions are: i) The pure integration is the simplest method but it has some problems like offset, initial rotor 
position, poor performance at low speeds, and stator resistance variation, ii) The LPF estimator was intended 
to solve pure integration problems but its dynamic response is bad and has drift in the flux magnitude, iii) The 
modified Integrator has good steady state and transient responses better than LPF and pure integration methods 
but it is a complex method and requires tuning for both LPFs and PI controller, and iv) The tunable LPF 
estimator has simple structure with excellent steady state and transient responses better than the other methods. 
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